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Growing complex metallic crystals, supported high index facet nanocrystal composites and tunable porosity 
metals, and exploiting factors that influence shape and morphology is crucial in many exciting 
developments in chemistry, catalysis, biotechnology and nanoscience. Assembly, organization and ordered 
crystallization of nanostructures into complex shapes requires understanding of the building blocks and 
their association, and this relationship can define the many physical properties of crystals and their 
assemblies. Understanding crystal evolution pathways is required for controlled deposition onto surfaces. 
Here, complex metallic crystals on the nano- and microscale, carbon supported nanoparticles, and spinodal 
porous noble metals with defined inter-feature distances in 3D, are accomplished in the solid-state for Au, 
Ag, Pd, and Re. Bottom-up growth and positioning is possible through competitive coarsening of mobile 
nanoparticles and their site-specific crystallization in a nucleation-dewetted matrix. Shape evolution, 
density and growth mechanism of complex metallic crystals and porous metals can be imaged during 
growth. 

Metallic nanocrystals with complex polyhedral shapes 13 have seen a surge in development 35 , owing to the 
lure of exotic and genuinely useful properties ranging from chemistry 3 to optics, and catalysis 6,7 to 
biomedical labeling, and high throughput combinatorial biotechnological 8,9 and other useful applica- 
tions where the nanoscale feature morphology and topology greatly influences performance 10 . The exquisite 
control now possible using solution-based methods 11,12 allows a myriad of complex shapes 13 using well-developed 
design rules 14,15 to control selective deposition onto specific crystallographic faces of the crystal using surfactants 
or ligands 11,14 . Crystals with complex shapes and composition have been made possible through the development 
of highly controllable protocols and design rules for their synthesis 5 in solutions. Standard methods require 
controlled seeding 2 and growth tuned by surfactants, with a strong dependence on the seed, which is sometimes a 
different material 4 . The dependence on solution-based methods limits the possibility of achieving similar and 
indeed new structures, some thermodynamically impossible in solution, in the solid state and on larger length 
scales. Solution-free, 'dry' syntheses are a necessity for some investigations and applications. Solution approaches 
offer bespoke shape design by relying on chemical engineering of homo- and heterogeneous phases or seeds, 
which must be single crystal 16 . Complex Platonic shapes 17 are often best controlled if the nucleation seeds 
themselves consist of significant complexity in shape and growth pathways were very recently observed when 
nanocrystal growth in solutions was monitored in real-time 4 . Understanding how nanoscale particles can assoc- 
iate during a complex decomposition process of a solid state precursor is a fundamental requirement to under- 
standing how the crystals grow 18 , and how the composition of the material can drastically affect the resulting 
morphology 19 , the ability to localise their deposition, or indeed write or pattern them onto surfaces. 

Since NPs and other nanoscale objects introduced into polymers can demix and coarsen into clusters through 
collective nucleation processes such as spinodal clustering 20 22 , they also prevent polymer dewetting from a 
surface 23,24 . These two factors inspired the hypothesis that shape controlled crystals from nanoscale to micron 
sizes could be controllably grown, positioned and the growth mechanisms monitored, directly on surfaces. 
Defining the role of nanoscale crystals and their influence on the crystallization of their coarsening or assembly, 
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is a major challenge for materials science 25 . The surface-growth cap- 
ability allows for growth pathways to be tracked in real-time using 
standard surface microscopies amenable to many laboratories. The 
formation of patterned nanoscale metallic structures by solid state 
dewetting 26 has proven useful for array formation and in porosity 
development 27,28 , but is often limited thin metal films in the first 
instance, where complex shaped crystals are not possible, and the 
metal content and resulting composite structure cannot be tuned as it 
can be with metal-polymer mixtures. 

Results 

We designed cyclic polymer and trimer-based complexes incorp- 
orating Au, Ag, Pd, and Re at the molecular level to allow pyrolytic 
demixing of physical mixtures of an organometallic derivative or 
metal salt, and polyspyro-, and spyrocyclic polyphosphazene (see 
Supplementary material). [N 3 P 3 (0 2 C 12 H 8 )3] or [NP(0 2 C 12 H 8 ) 2 ] n 
(0 2 C 12 H 8 = 2,2'-dioxy-l,l'-biphenyl) (Supplementary materials, 
Fig. SI) were synthesized and prepared that have X 5 -dioxydiaryl 
groups in the repeating units (polyspyrophosphazenes). The metal - 
lopolymer forms a bottom-up, spin-cast precursor; the metal is coor- 
dinated at a predefined position within the macromolecule. The 
pyrolytic approach and chemical decomposition processes these 
polymers undergo are summarized in Supplementary materials, 
Figs. Sl-3. 

For Au, the overall process leading to complex crystal shapes 
through demixing is summarized in Fig. 1. The cleaved metal atoms 
from early stage decomposition form NPs within the organic matrix 
(see Supplementary material, Fig. S4). Initial nucleation dewetting of 
the composite metallopolymer (Fig. 1A1), leads to the beginning of 
pinhole void formation (Fig. 1A2), whose complete beading is sup- 
pressed by the inclusion of metallic NPs 29 , observed by STEM to be 
distributed within a labyrinthine (spinodal-like) network of the car- 
bonizing polymer with a dominant length scale and interconnectivity 



between pinholes (Figs. 1B1, 1B2 and details in Supplementary 
material, Figs. S4,5). The pinhole density becomes constant after 
solidification of the carbonized polymer matter. 

Heterogeneous nucleation in the early stages (with spinodal char- 
acter, Fig. 1A2) occurs where the metallopolymer is liquid for a brief 
time above T g , allows mobile NPs to migrate solely toward the pinned 
rims of the pinholes, eventually ripening to a platelet seed (Fig. 1B3). 
For all processes in this work, the remarkable finding is that the 
pinholes form miniature reactors where crystallization and growth 
uniquely occurs. The process of spatially pinned seed formation and 
incubation is trackable using AFM, where in Fig. 1C, we observe the 
clustering 20 of NPs and the formation of 2D platelet seeds within the 
pinholes. We find that the subsequent crystallization of Platonic and 
complex shaped crystals uniquely from these pinholes (Figs. 1 A3 and 
ID), including geometrically perfect tetrahedra, decahedra and even 
single crystal metallic truncated icosahedra. A similar pathway lead- 
ing to specific crystal shapes is also demonstrated for Ag, Pd and Re 
(see Supplementary material, Fig. S6) using appropriately design 
macromolecular structures mixed with the corresponding metal salts. 

Figure IE demonstrates that halting of the pyrolysis indicates that 
the pinhole density is always greater than the density of crystals that 
eventual grow within them, until the vast majority of pinhole growth 
sites contains crystals. In this case, sufficient Au is present within the 
matrix to fill all available pinhole sites that form as a tessellation in 
the solidified matrix; the percentage of crystal-filled pinholes 
increases from 72% (at the stage where solidification is complete at 
~ 500°, to < 93% after full pyrolysis and annealing. Imaging after 
oxygen plasma etch removal of remnant organic matter confirms 
that the undulating features are not solely due to clusters forming 
under the polymer but from the instability throughout the matrix. 
The surface buckling is not due to spatial variations in surface tension 
during the brief period the metallopolymer is liquefied and impor- 
tantly, the films never dewet the substrate. Some definitive features 





Figure 1 | Tracked mechanism of solid state NP and complex crystal growth using Au. (A): (1-3) Pinhole formation in a solidifying metallopolymer 
during heating, forming a fixed array of pinholes leading to predefined density of shaped Au crystals. (B): ( 1 ) HAADF TEM images of gold NPs formed in 
the early decomposition stages. The lighter background features are the polymer matrix, bright points are the NPs. (2-3) Aggregating NPs migrate in- 
plane and from within the matrix to the pinhole acting as a nucleation point. Seeds form from Au NPs uniquely within pinholes. ( C) : AFM images of early 
stage growth where the nucleation of a crystal through NP coarsening within pinholes in the polymer is observed. While aggregation occurs at the 
nucleation pinhole, a single 2D crystal is first formed in all cases. (D): geometrically perfect micron-sized decahedron, tetrahedron and truncated 
icosahedron of single crystal gold grown in the solid state. (E): series of SEM images at frozen stages of crystal growth of Au crystals. The spatial density 
follows the pattern of pinholes on the surface once the seed is nucleated within a pinhole. Empty growth sites are shown in the outlined regions. 
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are observed during the crystal growth evolution. The nearest neigh- 
bor distance is typically related to the density of pinholes, and the 
density of pinholes defines the number of incubated crystals. This 
observation occurs with sufficient gold content in the macromole- 
cular precursors to fee the crystal growth. As will be explained below, 
crystal growth depletes the region around the crystal of Au NPs, and 
thus limits the development (density) and size of next-nearest neigh- 
bor crystals and reduces the density, leaving empty pinholes 
(Fig. IE). This phenomenon is also observed for Pd and Ag NP 
and crystal growth systems reported here. 

These crystals ripen from the pinholes uniquely, and EDX, XRD 
and XPS for these systems (see Supplementary material, Figs. S7-9) 
confirm the phase-pure, single crystal metallic composition of each 
metal within the matrix without surface bound contaminants of 
growth-influencing organics or carbon residues from the pyrolysis, 
as seen in HRTEM analysis of high index facet surfaces (see 
Supplementary material, Fig. S10). The resulting crystals formed 
are free form a carbon coating and with a defined areal density. 

Detailed AFM examination shows a consistent link (for a given 
metahorganic ratio) between the thickness of the deposit and the 
density of pinholes (as expected) but also the shape of the resulting 
crystals; thinnest deposits (60-150 nm) give tetrahedral growth (see 
Fig. 2), thickest deposits (> 3 um) resulted in pathways to metallic 
truncated icosahedral crystal growth. 

As we show in Fig. S7 and Fig. 2, cuboctahedra and tetrahedra are 
tracked at various stages of 3D growth. The crystal can be grown on 
almost any surface, or any shaped surface and the crystal shape 
evolution is maintained from nanoscale NPs up to micron scale 
crystals. While solution based synthesis offer mass production for 
a given size dispersion 30 , it approach requires additional methods 
often utilizing surface charge, potential gradients, or other deposition 
methods such as direct writing, droplet moulding 31 , and parallel 



printing to control position and location of individual or assembled 
materials. Pinning the entire crystal growth at a defined location also 
allows the growth method to be more amenably followed. 

We investigated the growth evolution of the crystal in detail (in 
terms of crystal growth and symmetry rather than ripening or coar- 
sening kinetic of the source NPs) and several pathways are outlined 
in Fig. 2 for Au; the types of crystals found for Re, Pd and Ag are 
shown in Supplementary material, Fig. S6 for comparison. The 
ripening steps of the tetrahedral Au microcrystals can be monitored 
from nanoscale clipped triangular platelets 32 (initial growth) to full 
tetrahedra (final structure), shown in Fig. 2A(l-7). As in solution 
systems, the tetrahedra evolve asymmetrically, where the edge length 
is first defined by forming a vertex-clipped 2D platelet delineated by 
dominant {111} 4 facets. Here, this shape maintains the same geo- 
metry from nm sizes through to micron scale as a 2D platelet due to 
being pinned at its incubation pinhole (NP nucleation and crystal- 
lization) site. The platelet then evolves to three-fold rotational sym- 
metry to a triangle with fully formed vertices, but always maintains 
its initial 2D footprint area (even when, by this method, the incuba- 
tion pinhole is larger than neighboring pinholes). Subsequent growth 
after the vertices are defined, upfills the structure until a geometric- 
ally perfect tetrahedron is formed 212 . The net effect contributes to the 
growth of clipped 2D platelets, observed many times in solution 
based approaches from nm to micron scale before forming the tri- 
angular platelet for continued growth (Fig. 1 and S7). 

On the nanoscale, penta-twinned seed NPs are found to grow to 
nanoscale icosahedra, but subsequently evolve to micro-scale deca- 
hedra (pathway Fig. 2B) as was recently suggested using heterome- 
tallic plasmonic labelling experiments on the nanoscale 4 . This 
observation confirms our hypothesis that the crystals fed from lat- 
erally mobile NPs evolve until all low energy {111} facets are dom- 
inant. Importantly, the method of tracking ripening growth allows a 
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Figure 2 | Tracking solid state crystal growth. (Al-7) The evolution of clipped triangular platelets and vertex formation to form a complete triangular 
platelet; subsequent growth results in a geometrically precise tetrahedron. (B) Penta-twinned NP that evolves to 20-fold icosahedron on the nanoscale. 
(Bl-3) The evolution of icosahedra to micron scale decahedra by anisotropic growth of clipped decahedra with subsequent outfilling from a pinned 
substrate position. (Cl-4) Evolution of tetrahedral on the micron scale to truncated bitetrahedra and eventually (via schematic on the left) to truncated 
icosahedral crystal form. Unmarked scale bars = 1 um. 
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true assessment of the shape evolution from the feedstock of metal 
atoms and nanoparticles that move towards the growth site in the 
dewetted, solidified polymer, without the complication of attach- 
ment of separately formed shapes. For example, oriented attachment 
of two tetrahedra to form a bitetrahedra would incorrectly describe 
the shape evolution process, which is difficult to monitor and control 
in solution syntheses for a single growing crystal. 

In contrast to recent findings 4 however, the penta-twinned NPs 
during crystallization are not found to evolve to tetrahedra (pathway 
B to C, Fig. 2). We consistently find that tetrahedra are observed to 
grow via the addition of extra twin planes by a different pathway - the 
primary difference compared to solid state thin film dewetting or all 
solution syntheses is the pinning of the crystallization site. 
Decahedra are not observed to evolve to 20-fold twinned icosahedra; 
tetrahedra on the micron scale, by continued preferential ripening 
from NPs within the spinodally decomposing metallopolymer, 
evolve by the addition of pairs of new triangles in place of the original 
six edges of the tetrahedra 1333 . The structure then evolves to non- 
uniform truncated icosahedra and then to perfect truncated 
icosahedral fullerene shapes (pathway C1-C4), a pathway not prev- 
iously seen for any size of metallic or other inorganic crystal, and its 
growth was tracked from a 2D nanoscale platelet seed (Fig. 2A1-7 + 
Cl-4) in the solid state. In all cases where the metastable metallopo- 
lymer thickness results in specific shapes (Fig. 1), micron sized 
shapes grow from specific nanoscale shapes, beyond thermodynamic 
limitations of solution growth. 

Tracking the growth mechanism allows and shape evolution path- 
ways offers the potential to create a variety of new structures for 
materials that can phase separate in this manner, including a range 
of organometallic precursors, and possibly others unobtainable by 
standard solution-based methods. We investigated the possibility of 
using nucleation and phase demixing/spinodal clustering, to create 
free-standing and spin-cast film composites for catalysis. Demixing 
and spinodal-type processes are dependent on many factors, such as 
the NP loading, thickness, m.w. and viscosity etc. Due to chemical 
changes in the polymer during pyrolysis in air, i.e. the carbon support 
(see Supplementary material, Figs. S2,3), its solidification, and the 
ripening of NPs (larger NP sizes with greater interparticle spacing), 
these factors dynamically change. Scaffolds for nanoparticles have 
been previously detailed for colloidal NP systems, and protocols for 
important dispersion of nano-objects in supports are well estab- 
lished. We find that out method allows for a definable switch in 
morphology. When the mixture ratio is 1 : 10 (larger polymer frac- 
tion), we find that mildly faceted NPs follow the surface pattern of the 
pinholes, but on a thick carbon residue (Fig. 3A1). For thicker even- 
tual deposits, the NPs are found to be smaller in size but do not 
develop faceting and are spherical (Fig. 3A2); for the thickest (up 
to 30 um, near-bulk) deposits we find NP-embedded carbon sup- 
ports for all metals examined. For Re (Fig. 3B), high-index faceted 
NPs (shown in more detail in Supplementary material, Fig. S10) form 
within a carbon support, where it characteristically forms via a 
layered growth process as the crystal emerges from the pinhole. Re 
crystals are unique in that they do not form defined polyhedral 
shapes on the micron scale. In the Supplementary material, Fig. S6, 
the larger crystal of Re form a layered structure, where the layers 
clearly form in a nautilus-like fashion from crystallization of coales- 
cing NPs with lower mobility. The crystals rotate as they grow from 
the pinholes giving the characteristic high surface area, layered mor- 
phology. Analagous to the uniqueness of the Re crystallization by 
comparison to Ag, Au and Pd, the corresponding early stage NPs 
shown in Fig. 3C, are high index faceted NPs with a narrow size 
distribution centered at a typical NP size of ~ 3 nm (Fig. 3D). XPS 
measurements confirm that the metallic NPs are freely accessible in 
the porous support and are not coated with carbon or other con- 
taminants. The method allows direct infiltration of NPs with a car- 
bon host, which is made porous through the formation of pinhole 
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Figure 3 | Tuning crystal size and forming NP-loaded layered graphitic 
composites. (A) SEM images of pseudo-spherical, partially faceted Au 
crystals formed from a (1) 3:1 trimer:gold mixture (2) 5 : 1 mixture ratio 
and (3) 10 : 1 NP(0 2 C 12 H 8 ) 3 /AuCl(PPh 3 ) ratio; the inset EDX map 
confirms phase pure Au. As the carbon content increases, so does the 
deposit thickness and the propensity for higher internal crystal density, 
with decreasing NP size (less Au). (B) Direct formation of high-index facet 
Re NPs within a carbon support. (C) Size dispersion of high index Re NPs 
within the graphitic matrix. (D) Raman scattering of the surface carbon 
deposit confirming conversion to graphite where D and G bands and 
harmonics related to crystalline defective graphite are tracked in time. 
(E) A ~ 1 mm thick deposit on Si0 2 , functionalized using Rhodamine (B) 
dye to show the morphology of the NP-loaded film. 



incubation sites for the growing metallic crystals. The formation an 
N0 2 , CO and C0 2 gases (see Supplementary material, Fig. S2 and 
S3), is also likely to contribute to the internal porosity of the carbon 
support. 

Raman scattering investigations (Fig. 3D) of the carbonization at 
several stages shows the emergence of convoluted scattered phonon 
processes as definitive D and G-bands, describing a mixture of 2D 
and 3D graphite eventually forming of a layered 3D graphitic foam 34 
(analysed in Supplementary material, Fig. SI 1) which is specifically 
due to the solid carbonaceous cyclomatrix formation (see Supple- 
mentary material, Fig. S12). The spectra evolve during pyrolysis to 
fingerprints of layered graphite, with remarkable similarity to layered 
graphene flakes with 'defective' curvature 35 . Using AFM, we find that 
the thickness significantly increases with higher carbon content, 
forming a graphitic 3D foam, from ~ 30 urn to over 1 mm in thick- 
ness on a surface, loaded with NPs (Fig. 3A3). In this case, subject to a 
specific further investigation, the overarching mechanism allows for 
a defined NP loading density based on the structure of the metallo- 
polymer and its phase demixing. This provides a tunable, single-step 
alternative to NP infiltration techniques and avoids agglomeration 
while maximizing spatial dispersion due to the structure and decom- 
position characteristics of the metallopolymer. During thermal 
decomposition, the mixture ratios that provide such composite 
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Figure 4 | 2D and 3D porous metal film crystallization. 

(A-C) Hierarchically porous Au formed by a spinodal decomposition of a 
[NP(0 2 C 12 H 8 )] n /AuCl(PPh 3 ) mixture with high Au content. HRTEM 
(inset) confirms the mesoporous metal structure. The smallest instability 
or feature wavelength occurs for the thickest films. (D), (E) Porous Ag and 
Pd thin films formed without dealloying, confirmed by back- scattered 
secondary electron imaging (inset) and (F), thick porous Au layer with 
consistent 3D inter-feature separation, cast on oxidized Si wafer. The insets 
are radially average FFTs of the porous metals confirming constant feature- 
feature separation in all directions. 

structure can have the size and density of NPs tuned. Short times 
provide dense, small NPs and longer heating times provid fewer, 
larger crystals. Figure 3E shows a how a thin film porous carbon film 
that is interspersed with metallic crystals can be uniformly dye- 
functionalized. 

Finally, we show that precursors with high metal content allow a 
facile one step route to pure porous metals and foams with a spinodal, 
3D disordered gyroidal morphology (Fig. 4), either as thick, 3D 
porous films, or as 2D porous metal coatings. Porous metals are an 
important class of materials used in ultra-sensitive SERS substrates 36 
and have unique properties such as high stiffness, low density, gas 
permeability and thermal conductivity for catalysis, hydrogen stor- 
age and acoustic insulation. Preparative methods for metal foams to 
date most commonly involve de-alloying of bimetallic Au x Ag„ 36,37 . 
Recent methods have achieved some control in porous metal forma- 
tion by spinodal and melt-induced dewetting of metallic thin films 27 . 
Controlling the distance between metallic features in stochastically 
dealloyed foams has been theoretically examined, but experimental 
verification of dealloying is limited for other metals; tuning the 
length scale between high aspect ratio features (proximity of sharp 
edges) is paramount for the enormous application possibilities of 



these materials. In SERS, the control of inter-feature width and 
curvature (edge definition) is crucial for maximal EM-field enhance- 
ment 38 . Maintaining a given porosity and porous structure from 2D 
to 3D, has not been achieved for Au or other metals, even though a 
very wide variety of metals can be dewetted in nanoscale features 
through localized heating 27 or ion bombardment 39 . This capability is 
useful in cases where dealloying of phases other than Au x Ag y is 
desirable, but not possible without a good contrast in standard reduc- 
tion potentials for the alloy components. By changing the metakpo- 
lymer ratio to a high metal content (see Supplementary material), the 
decomposition mechanism allows the spinodal-like patterning to 
occur for all metals examined, as confirmed through high resolution 
electron and surface probe microscopies in Fig. 4. The bottom up 
growth of the porous structure method differs from heating of thin 
metallic films, which is a top down dewetting process which inherent 
limitations in the degree of porosity that is possible. The formation of 
the characteristic percolating lateral topography with a 3D labyrinth- 
ine overall morphology 29 , initiated at the substrate interface, then 
maintains a characteristic and constant feature-feature distance in 
the metallic mesoporous morphology in 3D, even for layers that are 
several microns thick (Fig. 4F). In Fig. 4(A-E), the rings in Fourier 
space do not strictly refer to the formation mechanism, although they 
can often be used as signatures of spinodal decomposition in systems 
where this has been established; here, they confirm the random ori- 
entation and the degree of consistency of the inter-feature distances 
in each case, for both 2D films and thick, 3D porous layers, without 
adverse effects on SERS from remnant Ag often present in dealloyed 
materials 40 . 

Discussion 

Consistent with nucleation and demixing, all the remaining polymer 
between the pinholes in the thermalized deposits on surfaces 
becomes depleted of NPs as the lateral coarsening-feeding of the 
crystal growth continues; crystal egress continues 'layer-by-layer' 
from ripening NPs within the polymer structure, and also from 
NPs located on and near the surface. The younger ripening seed 
crystals form with good monodispersity prior to shape evolution, 
and their density is commensurate with metallopolymer thickness 
due to the uniform metallic content in the precursors, a denser dis- 
persion of NPs (of similar size and shape) is characteristically found 
from thicker deposits (see Supplementary material, Fig. S13) 28 . The 
mechanism of pinhole formation is progressive and as the carbona- 
ceous materials solidifies and hardens the surrounding materials, the 
resulting crystals have a variable size and degree of development. The 
preference for Au NP clustering indicates a strong (preferential) 
interparticle interaction 41 , where small seeds coalesce with a spatially 
pinned growing seed (in the pinholes) driven by phase separation 
rather than just coalescence 20 . This coarsening in the solid state is not 
likely to be dominantly controlled by thin-film 2D confinement as 
was found for aggregation of inclusions in simple PS thin films. The 
NPs also diffuse faster than aggregates of NPs tethered by polymer 
chains 42 , which do not ripen to a single geometrical crystal. 

It should be noted that the final diameters of the pinholes (crystal 
growth 'reactor' sites) remain unchanged after carbonization solidi- 
fication. The pinholes as nucleation points become pinned, essen- 
tially defining where crystal growth will only occur. In related solid 
state approaches such as dewetting of thin metallic films, break up to 
form NPs or meandering porous networks does not typically allow a 
reforming of a defined geometrical crystal shapes. Here, metal atoms 
migrate by spinodal clustering to the seeding crystal in a manner that 
preserves perfect crystal symmetry without changing the geometry 
during ripening and growth. Analysis of NP density from STEM 
images confirms a dominant 2D phase separation process (see 
Supplementary material, Figs. S4 and S13); for thicker films, NPs 
can migrate from within the decomposing mixture to the nucleation 
site (pinholes open the deposit from air to substrate). Each final 
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shape of crystal maintains a near-equal volume and edge length (see 
Supplementary material, Fig. S7) in spite of nucleation and growth 
being progressive rather than instantaneous and coupled. 
Consistently, all equivalent shapes finish with similar outer dimen- 
sions and overall volume. 

In solution, large jumps in crystal size are typically accompanied 
by evolution to thermodynamically lower energy shapes, such as 
spheres. The energy minimization has been demonstrated experi- 
mentally and theoretically 13 for nanoscale systems in solution. 
These calculations predicted a size threshold where the energy bar- 
rier for structural change becomes too great. The solid state decom- 
position of a molecular level building block metallopolymer 
undergoing phase demixing allows each crystal shape to evolve 
and maintain its complex geometry over a range of length scales - 
the shape evolution pathway is trackable. The density of pinholes is 
defined by the mixture ratio and molecular weights, and this ulti- 
mately defines the spatial density for a given crystal shape, for any 
stage of the shape evolution, which can be halted by quench cooling 
on a heat sink. Figure 1 showed stages of shape evolution and spatial 
density for triangular platelets that eventually grow to tetrahedra 
(Fig. 2), where the pinhole density is defined prior to the final crystal 
density (for sufficient metal content). 

The porous metals form by the combination of a spinodal-type 
instability but dominantly, a crystallization of a high NP coarsening, 
mimicking a spinodalfy decomposed polymer with metal NP inclu- 
sions 43 . The 'inverted' porous metal formation occurs in N 3 P 3 (0 2 
C 12 H 8 ) 3 /AuCl(PPh 3 ) and [NP(0 2 C 12 H 8 )] n /AuCl(PPh 3 ) mixtures, 
because a strong Au-arene interaction (between the AuCl(PPh 3 ) 
molecules and the biphenyl rings) allows regular ordering of the 
Au centers along the NP(0 2 C 12 H 8 ) groups of the polymer chains, 
or around the cyclic molecules of trimers. The NP association length 
however, is difficult to follow inside a carbon support during forma- 
tion. The instability development and porous metal topology could 
allow a wide range of functional metallic and oxide porous materials. 
The definition of the porosity and equidistance metal-metal dis- 
tances in porous metals is linked definitively to the molecular struc- 
ture of the precursor and the demixing/decomposition mechanism is 
likely to be adaptable to many metal-containing inorganic and metal- 
lic precursors. Defined deposition and porosity in porous Pd for 
example, as a pure metal or as a NP-containing composite augers 
well for catalytic application used in Li-0 2 batteries 44 and Suzuki, 
Heck and Negishi reactions, and gas sensing as pertinent examples. 
The dynamics of high metal loading and crystallization might also 
provide opportunities for periodic and complex, ordered porous 
metallic structure formation. 

Fundamentally, the observation of competitive demixing/dewet- 
ting of the polymer phase during decomposition the kinetics of crys- 
tallization from mobile NPs into complex shapes crystals and metals 
with tunable porosity will need to be examined in further detail. 
Some features have been established: the mobile NPs form specific 
crystals (when the minimal gold content and thickness of the deposit 
are defined) while the polymeric phase is liquid, since negligible Au 
content is found in the surrounding polymer phase after carboniza- 
tion (except for the mixture ratios that form NP-loaded carbon 
composites). NP association on and within polymers has been tested 
using Avrami and related kinetics, while nucleation and spinodal 
decomposition mechanisms (material independent) can accurately 
describe the pinhole formation and labyrinthine structure formation 
in simpler systems. Future work might identify the correct mech- 
anism involved in the competitive scenario shown here for greater 
composition and structural control of the resulting materials and 
patterns, particularly with progressively increasing viscosity (which 
occurs with an increase in temperature, but due to solidification) of 
the material during crystal growth. 

This bottom-up single-step approach for molecular level metallo- 
polymer building blocks that crystallize into a variety of defined 



structures and shapes provides a single step growth method without 
ligands or surfactants in solution to influence growth or shape. The 
association of NPs during phase demixing and the decomposition of 
the polymeric phases, which can be monitored, ultimately define the 
resulting structures into complex polyhedral crystals, metal-loaded 
porous carbon supports or porous metals. The method allows almost 
all important polymorphs of metallic crystals to form including an 
ability to monitor and delineate the growth process. Additionally, the 
growth pathways can be observed at all stages, where lateral demixing 
and subsequent ripening process give a well-defined crystal through 
a trackable shape evolution pathway. The areal density of the crystals 
depends on the thickness of the deposit, the structure of the macro- 
molecules, the metal loading and the competitive interplay between 
physical and chemical decomposition pathways allowing a transition 
between surface-bound individual crystals, carbon supported high 
index NPs, and tunable porosity metals 45 . This is made possible by 
demixing of a metal-containing macromolecular complex and co- 
processes of NP seed formation and incubated crystal growth in a 
graphitic matrix. When other materials are chosen, the approach has 
potential for multifunctional materials, liquid and water-free inor- 
ganic crystal synthesis for charge storage applications, plasmonics 
and SERS involving large scale porous metal films on any type of 
surface, photovoltaics 46 , and sensing when transition metals are used 
to form patterned oxides that can be deposited or even written using 
AFM onto surfaces. Scope also exists for supracrystallization of vari- 
ous crystal sub-units into complex functional assemblies which can 
be readily observed during growth. For example, porous loaded car- 
bons improved by further development may be applied to catalytic 
composites using small NP loaded carbons for conversion processes 
in catalysis and energy storage conversion processes where the 
placement NPs into support matrices without agglomeration and 
shape changes is possible during synthesis. This approach can be 
readily extended to a very wide range of transition metal oxide, 
metalloid and semiconductor crystal growth, composite formation 
and porous assemblies. 

Methods 

Solvents used were dried a purified using standard procedures. The starting materials 
[NP(0 2 C 12 H s )] n (with a m.w. of 600,000) and [N 3 P 3 (0 2 C 12 H 8 ) 3 ] (1) were prepared 
using the reported methods. The mixture complexes were prepared in high yield 
(96%) and evaporating the solvent. The purity of the product was verified by 31 P- 
NMR. Details on the mixtures, and their composition are as follows: 

Preparation of the mixtures. AuCl(PPh 3 )/[NP(0 2 C 12 H s )]„. The AuCl(PPh 3 ) 
complex was prepared in high yield by adding AuCl (Aldrich) to a dichloromethane 
solution of PPh 3 (1:1 molar ratio) and evaporating the solvent. Table 1 summarises 
the experimental conditions for the preparation of the mixtures. 

AgPPh 3 [CP 3 S03]/[NP(0 2 C 12 H s )] 3 . AgPPh 3 [CF 3 S0 3 ] (0.3795 g, 9.5 X 10" 4 mmol) 
and [NP(0 2 Ci 2 H s )] 3 (0.5130 g, 7.46 X 10~ 4 mmol) for a molar relationship 1 : 1; 
AgPPh 3 [CF 3 S0 3 ] (0.074 g, 1.86 X 10-" mmol) and [NP(0 2 C 12 H 8 )] 3 (0.4234 g, 6.16 
X 10~ 4 mmol) for a molar relationship 1 : 5 and AgPPh 3 [CF 3 S0 3 ] (0.0594 g, 1.5 X 
10~ 3 mmol) and [NP(0 2 Ci 2 H 8 )] 3 (0.4064 g, 5.9 X 10~ 4 mmol) for a molar rela- 
tionship 1 : 10, were stirred in a CH 2 C1 2 solution at room temperature for 24 h. After 
this the solvent was eliminated under reduced pressure and the solid dried under a 
vacuum at room temperature. 



Table 1 | Summary of the experimental conditions forth 
tion of the mixtures (n = 3 or polymer) 


eprepara- 


Mixture Polymer/Au 


Trimer/Au gofAu 


gofTrimer g 


of Polymer 


1 3:1° 

2 

3 

4 


0.1 1 
3:l c 0.004 
5:l c 0.0858 
10:1= 0.027 


0.017 

0.0238 

0.3865 


0.15 


"Film; 

b Power method. 
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PdCl 2 /[NP(0 2 C i2 H 8 ). The respective molar mixture of Palladium chloride {0.1 g, 
0.56 mmol for the 1.1 ratio mixture) and {0.14 g, 0.789 mmol for the 1 : 5 ratio 
mixture) and the trimer [NP(0 2 C 12 H 8 )] 3 {0.4 g, 0.58 mmol for the 1.1 ratio mixture 
and 1 .94 g, 2.82 mmol for the 1 : 5 ratio mixture) were dissolved in CH 2 C1 2 , stirred for 
24 hours and evaporated to dryness under vacuum at room temperature. 

K[Re0 4 ]/[NP(0 2 C 12 H s ). The mixtures of K[Re0 4 ] {0.1 g., 0.354 mmol) and the 
cyclotriphosphazene [NP(0 2 C 12 H 8 )] 3 {1.19 g., 1.73 mmol) for the 1 : 5 molar ratio 
and of K[Re0 4 ] (0.03 g., 0.103 mmol) and [NP(0 2 C 12 H 8 )] 3 (0.7 g., 1.018 mmol) for 
the 1 : 10 molar ratio, were dissolved in dichloromethane, stirred for 24 hours and 
evaporated under vacuum until dry. 

Powder. The respective molar mixtures above are each dissolved in THF, stirred for 
24 h and evaporated to complete dryness in a vacuum at room temperature. 

Deposition of the mixtures on Si or Si0 2 . A drop of CH 2 C1 2 solution containing the 
mixtures (solutions/ suspensions in the range 1 X 10" 3 -5 X 10~ 4 g/mL) was cast over 
either a silicon or oxidized silicon surface and the solvent evaporated. The Si or Si0 2 
wafer (with 400 nm of thermally grown oxide) was subsequently pyrolyzed (see 
Supplementary materials, Fig. SI) at 800 °C using a predefined temperature program 
(Daihan oven model FHP-12) where the polymer samples were heated at a rate of 
10°C min -1 from room temperature to 800 'C under a constant air flow of 120 mL 
min -1 and subsequently annealed for 2 h. The resulting depositions give a uniform 
areal coverage for each mixture. 

Analytical techniques. X-ray diffraction (XRD) was carried out at room temperature 
on a Siemens D-5000 diffractometer with 0-20 geometry. The XRD data was collected 
using Cu-Kot radiation {40 kV and 30 mA). Scanning electron microscopy (SEM) 
and energy dispersive X-ray analysis (EDX) were acquired with a JEOL 5410 SEM 
with a NORAN Instrument micro-probe transmission microscope and with a Hitachi 
SU-70 FESEM operating at 10 kV equipped with an Oxford Instruments X-max 
50 mm 2 solid-state EDX detector for elemental line scans and mapping. AFM 
measurements were acquired with a Veeco Instruments Enviroscope operating in 
intermittent contact mode. X-ray photoelectron spectroscopy (XPS) was acquired for 
deposits on Si and Si0 2 surfaces in normal emission geometry. A VG CLAM 2 
spectrometer with an Al Ka X-ray source (1486.6 eV), operating at 250 W, were used 
several centimeters from the samples. The spectra were normalized to the total 
electron yield to correct for small differences in sample positions and X-ray source 
intensities. The resulting spectra were fitted by using a Shirley-type background and 
symmetric Voigt functions with variable Gauss and Lorentz contributions. 

Image analysis of feature topology. AFM and optical microscopy were used for 
analysis of patterns, topologies and morphologies of the NP coarsening and overall 
crystallization process. For in-situ examination of crystal growth and cluster size 
formation, optical microscopy data and AFM data was acquired during pyrolysis up 
to the spinodal feature development stage. Both techniques were used to follow the 
variation of the in-plane morphology, i.e. the direction of NP migration and feeding 
of crystal growth and patterning. Height variations of patterns and crystal were 
obtained by AFM and compared to SEM. 

Structure factors used for pattern determination and crystal/pinhole relative 
positions on films were acquired by thresholding digital images and obtaining radially 
averaged fast Fourier transforms (FFTs). The maximum intensity of the structure 
factor S(k) from FFTs gives the wavelength of surface morphology in the case of 
spinodal-like pattern formation, where the wavelength X = 27t/k. Random distribu- 
tions of nanop articles and morphologies without any definitive wavelength in posi- 
tion within the morphology give a decaying intensity from the maximum intensity 
(000) central mode. Nanoparticle and crystal coverage within decomposed films from 
SEM and high resolution TEM images were determined using ImageJ analysis soft- 
ware for particle size and spatial distributions, accounting for shape anisotropy. NP 
cluster densities were measured using optical microscopy. 
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